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£"•»« Abstract: In this paper we show that the intuitive guess that the geometric scaling behaviour should be 
violated in the case of the running QCD coupling, turns out to be correct. The scattering amplitude of 
\Q the dipole with the size r depends on new dimensional scale: Aqcd, even at large values Y = ln(l/x) and 
V.Q / = In (as (r 2 ) /as fl/Q 2 ,)). However, in this region we found a new scaling behaviour: the amplitude is a 
^ function of £ = Y I. We state that only in the vicinity of the saturation scale Q s (as(Q 2 ) In (r 2 <5 2 ) < 1), 
T _— I the amplitude shows the geometric scaling behaviour. Based on these finding the geometric scaling behavior 
j> jthat has been seen experimentally, stems from either we have not probed the proton at HERA and the 
LHC deeply inside the saturation region or that there exists the mechanism of freezing of the QCD coupling 



$-h constant at r 2 rj 
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1. Introduction. 

The geometric scaling behaviour of the scattering amplitude gives an example of the prediction that based 
on the most fundamental features of the high parton density QCD. It means that the scattering amplitude 
of the dipole in the saturation region is a function of one dimensionless variable r = r 2 Q 2 (Y; b) instead 
of being a function of dipole size (r), energy (Y = ln(l/a;)) and impact parameter (6). Q s (x,b) is the 
new dimensional scale ( saturation momentum) which absorbed entire dependence on energy and impact 
parameter of the amplitude. The existence of this scale and its appearance from the non-linear evolution at 
low x is the most fundamental theoretical result of both the BFKL Pomeron calculus [1-5] and the Colour 
Glass Condensate (CGC) approach [6-9]. The idea of the geometric scaling behaviour is very simple: r is 
the only dimensionless variable in the dense system of partons. At the moment we have a general proof 
of the geometric scaling behaviour of the amplitude in the saturation region [10], two examples of the 
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analytically solved non- linear Balitsky-Kovchegov equation with simplified kernels [11,12] which show the 
geometric scaling behaviour and the proof that this behaviour is a general property of the linear dynamics 
in the vicinity of the saturation scale [13]. 

However, analyzing these arguments and proofs one can see that they are related to the case of fixed 
QCD coupling constant as far as the behaviour of the amplitude in the saturation domain is concerned. 
Indeed, for frozen QCD coupling we do not have other dimensional parameters but Q s . For running QCD 
coupling the situation is not so clear since it brings the second domensional scale: Aqcd, and the role of 
this scale has to be studied in the saturation region. 

The paper presents such study in the case of the simplified BFKL kernel (see Ref. [11]). We show 
that the geometic scaling behaviour is violated for the running QCD coupling and the amplitude does not 
depend on the one variable r. It turns out that the amplitude depends on a different variable 

AN ( gg { r2A QCD) \ 4N C 
C = —r- Yhi -, ^ where a s {r ) = ^- (1.1) 

b \*s(V Q cn/Qi)) b\n[l/[r^ QCD )) 

with b = lliV c /3 — 2Nf/3 for number of colours N c and the number of flavours Nf. One can see that £ 
depends on both dimemsional scales: Q s and Aqqd- 



2. General approach: behaviour of the scattering amplitude in the vicinity of the 
saturation scale for running QCD coupling 



The nonlinear Balitsky-Kovchegov equation for the scattering amplitude of the dipole with size r has the 
following form [8,9]: 

l)A{rA:h) f d2ri K ^ rur2 ) x l N ( ru Y;b-lf 2 )+N(r 2 ,Y;b-l-r 1 )- N (r, Y; b 



2 



dY / 2vr 



v j ry.Y-.b - lr 2 ) N (r 2 ,Y;b- ^ } (2.1) 



where Y = ln(l/x) is the rapidity of the incoming dipole; ./V is the imaginary part of the scattering 
amplitude and b is the impact parameter of this scattering process and f 2 = f ' — f\. The BFKL kernel 
K (r±,r 2 ) has the following form 



K(r;n,r 2 ) _ &s (r>) + » ( ^ - l) + \ ( - l) (2.2) 



This kernel takes into account the running QCD coupling and was derived in Re. [14]. In Eq. (2.1) as is 
the QCD coupling 

as[r} = i + ^mi n{R 2j^ = b My aqcd)) (2 ' 3) 
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and as = N c as/n- R is the arbitrary size (so called the renormalization point) which the physical 
observables do not depend on. 

In the vicinity of the saturation scale where r 2 ~ r\ ~ r| ~ 1 /Q 2 and we can consider that as (r 2 ) = 
( r i) = ols ( r i)- Indeed, choosing R = r we can see that 

as (r 2 ) ln(r 2 /r?) <Zn(r 2 A QCO ) 

1 + ^MrVr?) 



«» M) = . . ^li^l^A Qs M (2.4) 



In the vicinity of the saturation scale r 2 oc 1/Q 2 and condition | In (r 2 Q 2 ) | <C In (Q 2 / 'Aqcd) determines 
the kinematic region which we call vicinity of the saturation scale. Using this simplification the kernel of 
Eq. (^T|) looks as follows: 



K(r;ri,r 2 ) = a s (r 2 ) (2.5) 

The second simplification stems from the observation that for the equation for the saturation scale 
we do not need to know the precise form of non-linear term [2, 12, 15]. Therefore, to find this equation as 
well as behaviour of the amplitude in the vicinity of the saturation scale we need to solve the linear BFKL 
equation, but in the way which will be suitable for the solution of the non-linear equation with a general 
non-linear term. It is enough to use the semiclassical approximation for the amplitude N (r,Y;b), which 
has the form 

N A (Y,0 = e s ™ = e -(^ + (i-7(^)K + So ( 2 . 6 ) 

where £ = In (r 2 Q 2 (Y = Yq; &)) . In Eq. ( |2.6| ) we are searching for functions u (Y, £) and 7 (Y, £) which are 
smooth functions of both arguments in the following sense 

J Y (Y, « u> (Y, ; u/ 6 (Y, « u> (Y, ; 7 y (Y, « 7 (Y, ; 7 g (Y, £) « 7 (Y, £) ; (2.7) 

The BFKL equation near to the saturation scale looks as follows 
dN(r,Y;b) , 2N f d 2 n 



dY 



/Y2 
-^K(r;r 1 ,r 2 )x{N(r 1 ,Y;bj + N (r 2 ,Y;b) - n(t,Y;$)} (2. 



In Eq. J2.8D we assume that we are looking for the solution at b 3> r\ or /and r 2 . Substituting Eq. (|2i 
into Eq. ( |2.8| ) and taking into account that function (r 2 )-^ = exp(/£) is the eigenfunction of the BFKL 
equation, namely, 

a S {r 2 ) J ^-K(r-r u r 2 ) (r 2 / = a s (r 2 ) X (/) (r 2 ) f with x (/) = 2^(1) - ^(/) - ^(1 - /) 
where ifj(z) = dlnT(z)/dz and r(z) is Euler gamma function (2-9) 

we obtain that 

w(y,0 = a S (0 X(7(Y,0) (2-10) 

This solution has a form of wave-package and the critical line is the specific trajectory for this wave- 
package which coincides with the its front line. In other words, it is the trajectory on which the phase 
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velocity (v p h) for the wave-package is the same as the group velocity ( v gr ). The equation v gr = v p h has 
the folowing form for Eq. ( |2,10| ) 



v P h = a s (r 2 ) X ler) = -a s (r 2 ) x ilcr) = v gr (2.11) 

i 7cr 



with the solution j cr = 0.37. 



Eq. ( 2.11 ) can be translated into the following equation for the critical trajectory 



dj{Y) _ , c v X (lev) 10 s 

= = « s (6tt^ ( 2 - 12 ) 

with the solution 

^f<^y s £ = f>-*S (2.13) 

1 — 7cr 

where £ = In (Q 2 (Y = Y ; 6) /A| CD ) and £ s = In (Q 2 (Y &)/Q 2 (Y = Y , 6)) 

For finding the behaviour of the amplitude in the vicinity of the line given by Eq. (|2.13| ) one should 
expand function w(Y,£) and 7(Y;£) and find a deviation from the critical line of Eq. (|2.13[) . Replacing 
£ = £ s + A£ where A£ = In (r 2 Q 2 (Y, b)) and considering A£ ^ s one obtain 

N(Y,C) cc exp{( ^ ( ^ =6) Y + l- 7 ) Ag 
= «p/f%)y + 1 _ 7flp W 



e 



r 2 Q 2 (Y & )) |(1 - 7cr) (2.14) 

It should be mentioned that everything, except Eq. ( p,14| ), are not new and have been studied in 
details before (see for example Refs. [2, 12, 15]). We discuss them here for the complitness of presentation. 
Eq. (2.14) shows that the running QCD coupling leads to a different behaviour of the scattering amplitude 

-(l-7er) 



in the vicinity of the critical trajectory. Recall that for frozen as the amplitude N cc ( r Q s 
Concluding this section we would like to stress that we obtain the geometric scaling behaviour of the 
scattering amplitude to the right of the critical curve (r > 1) in the case of the running as- This result 
gives us a hope that inside the saturation region we can observe the geometric scaling behaviour as well. 



3. The non-linear equation with the simplified BFKL kernel for running as- 



The solution to the Balitsky-Kovchegov equation with the kernel of Eq. (2.2) has not been found. Following 
Ref. [11] we simplify the kernel by taking into account only log contributions. In other words, we would like 
to consider only leading twist contribution to the BFKL kernel, which contains all twists. Actually we have 
two types of the logarithmic contributions: In (r 2 IS^ CD J for r 2 <C l/Q 2 and In (r 2 Q 2 ) for r 2 > l/Q 2 . 
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3.1 r 2 < 1/Q 



In this kinematic region we can simplify K (r;r i,r 2 ) in Eq. (2.2) in the following way [11], since ri > r 
and r 2 = \r — r*\> r 

dr' 



d 2 r' K (r, n .r 2 ) -^7ra s (r 2 )r 2 / QCD — (3.1 



Introducing N (r, Y; b) = N (r, Y; b) / (as (r 2 ) r 2 ) we obtain 

dN (r, Y; b) f^Qon ( a s (r' 2 ) - a 2 {r' 2 ) 



dY 



J 2 *' 2 -^J- iV (r', F; 6) - ^-^N 2 (/, F; 6) \ (3.2) 



J" 2 QCD rfr' 2 s r \ 2 - I . As we have discussed 

in the previous section the form of non-linear corrections is not important here. One can see that the linear 
part of Eq. ( |3.2| ) gives the familiar GLAP equation in the double log approximation [16]. 



3.2 r 2 > 1/Q 2 

The main contribution in this kinematic region originates from the decay of the large size dipole into one 
small size dipole and one large size dipole. However, the size of the small dipole is still larger than 1/Q S . It 
turns out that as depends on the size of produced dipole if this size is the smallest one. It follows directly 
from Eq. (2^) in the kinematic regions: r ~ r 2 3> T\ ^ 1/Q 2 and r ~ r\ 3> r 2 3> 1/Q 2 (see Ref. [17] 
for additional arguments). This observation can be translated in the following form of the kernel 

Sr'K(ry) , f + « f (3.3) 

Jl/Q*<yjb) r i JVQ 2 s (Y,b) r 2 

One can see that this kernel leads to the (fi/Q2(Y,b) as ^) dri \ -contributions. Introducing a new 
function 

r 2 a?(r' 2 ) 

N(r,Y;b) = / dr' 2 \ 9 J N(r',Y;b) (3.4) 



one obtain the following equation 

dN (r, Y; b) 
dY 

Introducing a new variable 

- 2 - i J2 



N(r,Y;b) (l - N (r,Y;b)) (3.5) 



with £ = — In ^r 2 A 2 ^^ = — £ and new function (r, Y; 6) 

AT( r ,y ; &) = 1 - e -^ r ' Y V (3.7) 
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Y 



Saturation domain 



<P = % 



1 = Is 



Perturbative QCD region 



Figure 1: Kinematic regions for the dipolc 
scattering: £ = — In (r 2 Aq CD ) and £ s = 
In (Qs/^qcd) • The red line is the critical tra- 
jectory with the equation 32 ^ c Y = £ 2 (see 



Eq. (4.3)) on which 0(Y,£ = &) = fa. 



we obtain the following equation 



d 2 (j) (r, Y- b) 
dYdl 



(r,y ; i 



(3.8) 



4. Solutions to the simplified equation 



4.1 Initial and boundary conditions. 

The simplified BFKL kernel looks as follows [11] in uj and 7 representation (in double Mellin transform 
with respect to Y and £) 



( 1 



i for r 2 < 1/Q 2 ; 

for r 2 > 1/Q 2 ; 



(4.1) 



1-7 



Using this kernel for small values of r 2 ( r 2 < 1/Q 2 ) and the general formulae of Eq. ( 2.13j ) and 
Eq. ( |2.14 ) one can write the initial conditions at r = r 2 Q 2 = 1 . They are 



(y,l = 6;&) 



4^0 



(4.2) 



The critical line that gives us the energy dependence of the saturation scale has the form (see Eq. ( |2.13| ) 

32iV r 



Y 



F2 



(4.3) 



In Eq. ( p~3| ) we assume that £0 = 0. It means that we consider the scattering amplitude for the dipoles of 
all sizes smaller that tq = 1/Aqcd and the entire kinematic region can be divided in two parts: the region 
of perturbative QCD and the saturation domain (see Fig. Q). 
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4.2 Solution for > 1 

Searching for the solution to Eq. ([^) we start with finding the asymptotic berhaviour of 4> at large values 
of Y and I. We expect that (ft will be large in this region since that dipole amplitude tends to be close to 
unity due to unitarity constraints. Therefore, in this region Eq. fl3.8| ) degenerates to a very simple equation 

8Y81 K ' ' 

with obvious solution: 

0oo(Y,Z;&) = Yl + F(Y) + G (I) (4.5) 
where functions F and G should be found from the initial conditions of Eq. ( [4.2| ). 
The final solution has the form 

0oo (Y,l;b) = Y(l-l s ) - ~^> (e l - e l °) - J (e 21 - e 2l °) + <p (4.6) 
where/ s = ^Tn£ s . 



4 TU V J 2 



Therefore, we learned two lessons in this subsection: (1) the main problem with Eq. (3.8) is to satisfy 



the initial and boundary conditions; and (2) the asymptotic solution does not show a geometric scaling 



behaviour since even the simplest solution of Eq. (4.6) does not depend on the variable z = £ s — £. However, 



the solution of Eq. (|4.6j ) at £ — > £ s has the following form 

<MW) c ~ 6<<6 > h - l<h = <t> + l^oz (4.7) 

showing the geometric scaling behaviour. Hence, we can hope that the solution will show the geometric 
scaling behaviour in the vicinity of the saturation scale. 

On the other hand, the solution given by Eq. ( |4.6| ), leads to cf)^ < and, therefore, contradicts the 
unitarity constraints, leading to the negative imaginary part of the amplitude. Such behaviour stems from 
that terms in Eq. ( |4.6| ) which are responsible for the matching of the d(f)/dl at I — > 0. Hence we have to 
find a diffrent solution which has the same behaviour 0oo = Y I for ^ > 1. 

4.3 Traveling wave solution 



Eq. (|3.8| ) has general traveling wave solution (see Ref. [21] formula 3.5.3) which can be found noticing that 
(j>(Y,l;b) = <fi [r] = aY + bl;b) reduced the equation to 



di 

The general solution of Eq. (f4.§|) has the form 



l-ii) 
J <t>n 



ab^p- = l- e -«**> (4i 
dr] z 



r, = aY + bl (4.9) 
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where c,(po,a and b are arbitrary constants that should be found from the initial and boundary conditions. 



The initial conditions of Eq. (4.2) can be written in terms of Y and I variables as 



0(77 = oF + W,;6) = 0o ; <f>' n {v = aY + bl s ;b) = -j^ok (4.10) 



It should be mentioned that the variable rj is not the scaling variable z = In (r) = £ s — £ with £ s = y 32 ^ c V. 
One can see that we cannot satisfy the initial conditions of Eq. ( f4.10| ). Indeed, even to satisfy the first of 
Eq. ( 4.1C| ) we need to choose rj = on the critical line. As you see we cannot do this with a and b being 



constants. The second equation depends on Y, but not on rj, making impossible to satisfy this condition 
in the framework of traveling wave solution. 

If we try to find a solution which depends on z (0 (Y; r 2 ; b) = <ft (z; b)) we obtain the following equation 

1&N C 



(using the variable z = \ -^-^ z 



16iV c z d 2 <p(z;b) , dy\..u, , _-4>(z;b) 



f^Y dz 2 dz 2 



e -<p(m (4.11) 



Therefore, only in the vicinity of the critical line where \ 16 / Vc z <C y2Y we can expect the geometric 



b 

scaling behaviour of the scattering amplitude. It should be stressed that at large value of Y the region 
where we have the geometric scaling behaviour becomes rather large. Neglecting the first term in Eq. ( [4.11 ) 



we obtain the equation in the same form as for frozen as- It is easy to find the solution to this equation 



that satisfies the initial condition of Eq. (|4.2| ) . Actually, the condition y 16 ^ c z -C \/2Y can be rewritten 
as as (Qs) hi (r 2 Q 2 ^j <C 1 and it shows the region in which we can consider the running QCD coupling as 
being frozen at r 2 = l/Q 2 . 



4.4 Self-similar solution 



Generally speaking (see Ref. [21] formulae 3.4.1.1 and 3.5.2) Eq. (3.8) has a self similar (functional 
separable) solution <j)(Y,l;b) = (j)((;b) with ( see also Eq. ( |1.1[) ) 

C = Y{l-l s ) (4.12) 

For function (f> (£; b) we can reduce Eq. Q3.8]) to the ordinary differential equation 



(C - 2) + *M = 1 - e -«C*) (4.13) 

dQ z dQ 



The initial condition of Eq. (|4.2j ) can be rewritten in the form 

0(C = O;6) = O ; (4.14) 



d( 4™" 4 
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50 



Figure 2: Solution to Eq. fl4.13| ) for different values of ^ (C = 0; b) (from top to bottom ^ (( = 0; b) = 2 (black) ,0.2 
( blue), (red)). 



Generally speaking we cannot satisfy Eq. (glj) using the solution of Eq. d4.13|) since these conditions 
depend not only on £ but on extra variable Y. However, at large value of Y one can see that Eq. ( KW) 
degenerates to 

0(C = O;6) = ^; ^ (C = ° ;6) = 0; (4.15) 



which are consistent with the solution being the function of only In Fig. || we plot the numerical 
solution to the Eq. ( [4.13 ) at different values of 4>'^ (£ = 0; b). One can see that this solution is not sensitive 
to this value if it is small enough. It means that the £ scaling behaviour can start from rather small values 
of Y. Since the whole approach, based on leading log(l/x) contribution, can be trusted only at large values 
of Y we believe that £ scaling behaviour is a good approximation to the solution of Eq. (|3.8|) . 



4.5 Numerical solution 

It turns out that this believe was too optimistic. In Fig. || we plot the numerical solution to Eq. ( |3.S| ) with 
the initial condition given by Eq. (fi^). The main lesson that we can learn from these pictures is that the 
£- scaling behaviour can be reasonable approach but at unreasonably large values of I or/and Y. Indeed, 
at Fig. |^-b we can see that at large values of £ the solition depends on / only slowly. 

At large Y and I the exact solution is reasonable to compare with the solution of Eq. (4.13) (see Fig. |j). 
One can see the same pattern: they become close at large values of Y and I (£ and I). 



5. Conclusions 

In this paper we show that the intuitive guess that the running QCD coupling will violate the geomet- 
ric behaviour of the scattering amplitude, turns out to be correct. Indeed, we found out that the new 
dimensional scale: A.qcd, that brings the running as, enters to the amplitude behaviour even at very 
high energies. However, in the vicinity of the saturation scale (r 2 oc 1/Q 2 S ) we see the geometric scaling 
behaviour. This vicinity is determined by ocs(Q s ) m ( r2 Qs) ^ 1- I n other words, the geometric scaling 
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Fig. §b 



Fig. §a 

Figure 3: The exact solution of Eq. (3.8) with the initial conditions given by Eq. ( |4.2[ ) for function <f> {C,/l, I). Fig. ||-a 
gives the bahaviour of <fi at small values of £ and Z while Fig. ||-b shows the same behaviour in the region of large £ 
and Z. The value of <fo was taken <j>o = 0.1 
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Fig. §b 



Fig. g-a 

Figure 4: The exact solution of Eq. ( [3.8] ) with the initial conditions given by Eq. ( |4.2[ ) for function (Y, Z) versus 
^-scaling solution of Eq. ( 4.13 ) <// (£). In Fig. f|-a ln</> (C/Z, Z) and ln0' (£) are plotted at diffrent values of (. Fig. ||- 
b shows the same but for the amplitude N((/l,l) = 1- exp (-0 (C/Z, Z)) and A/ 7 (Y£) = 1 - exp (-</>' ((")). The 
calculations were performed at (j)Q = 0.1. 



behaviour of the amplitude remains until we can neglect the difference between as (r 2 ) and as (l/Q 2 ). 
In different way of saying, if we could find the mechanism that will freeze the running as on r 2 = 1 /Q 2 
the amplitude would show the geometric scaling behavior. However, in the framework of the leading twist 
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BFKL we did not find such a mechanism. 

For as(Q s ) hi (r 2 Q 2 ) > 1 the geometric scaling behaviour is violated and at very large Y and I = 
In ^ln (t 2 Aq(jj^J / In (Q^Aq CD ^ J the amplitude depends only on one variable ( (see Eq. (1.1)). 

From our point of view the fact that we see geometric scaling behavior experimetally, stems from either 
we have not probed the proton at HERA and the LHC deeply inside the saturation region or that there 
exists the mechanism of freezing of the coupling QCD constant at r 2 ~ 1/Q 2 . In practical calculations as 
is used to be frozen at some value of the momentum larger that l/r? r > Aq CD . In this case we would 

like to notice that if as(Qs) l n ( r /r^s) < 1 we still have the geometric scaling behaviour. For example 
in recent paper of Ref. [22] the value of r/ r is chosen from cts( r fr) = 0.7 or 1. For such value of r/ r we see 
that as(Q s ) In [r 2 r Q 2 A < 1 for Q 2 = 0.3 -f- 4GeV 2 covering the region of energy from RHIC to LHC. 
Therefore, in the CGC motivated model of Ref. [22] we do not expect to see any violation of the geometric 
scaling behaviour. 
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